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Abstract: The Energetic Gamma- Ray Experiment Telescope (EGRET) experiment onboard the Compton 
Gamma-Ray Observatory (CGRO) has provided the most precise measurements of the 7-ray sky to date. 
The EGRET measurements of the diffuse emission across the sky show an excess above 1 GeV. This "GeV 
excess" has been a topic of great debate and interest since its original discovery by Hunter et al. in 1997. 
We have modified the GLAST simulation and reconstruction software to model the EGRET instrument. 
This detailed modeling has allowed us to explore the parameters of the EGRET instrument, in both its 
beam-test configuration and in-orbit on CGRO, in greater detail than has previously been published. We 
have found that the GeV excess is significantly increased when previously unaccounted for instrumental 
effects are considered. We will present a new measurement of diffuse 7-ray emission in the inner Galaxy. 



Introduction 

The EGRET[1] telescope onboard the CGRO pro- 
vided the most detailed look at the 7-ray sky to 
date. Before launch, the EGRET collaboration per- 
formed extensive testing [2] to characterize the in- 
strument response to both charged particles and 7- 
rays. The instrument response functions found dur- 
ing these tests were the basis for constructing the 
EGRET exposure maps and thus the EGRET in- 
tensity sky maps. While there are references to a 
Monte Carlo simulation of the EGRET in a variety 
of EGRET publications [2, 3], it is clear from those 
publications that the simulations were preliminary 
in nature. 

The simulation environment developed for use 
with the Gamma-ray Large Area Space Telescope 
(GLAST) [4] provides a basis on which we have 
developed our simulations. The simulation frame- 
work is based on GEANT4 [5, 6] which is well 
tested. We have constructed a geometric model 
of the EGRET instrument based on the most de- 
tailed information available. We have also in- 
cluded parametric response models for the EGRET 
sub-detectors. Both the geometric and response 
models are incorporated into the framework. Fur- 
thermore, we have integrated the Burst and Tran- 



sient Source Experiment (BATSE) Mass Model of 
CGRO [7], used to estimate backgrounds in the 
BATSE instrument, with our model of the EGRET 
instrument. 

Characterization of the EGRET instrument was ex- 
tensive, consisting of a charged particle beam test 
and two 7-ray beam tests. However, it was not 
possible to test all configurations of the instru- 
ment, nor was it possible to examine the particle 
interactions outside of the detectors. The Monte 
Carlo framework we have constructed allows us to 
probe the instrument in detail and access informa- 
tion about the simulated events not available in the 
laboratory. Thus, we have been able to compare 
the differences between the EGRET in the beam 
test environments to the EGRET in its flight envi- 
ronment. Most notably, we have probed effects re- 
lating to the integration of the EGRET onto CGRO. 

The EGRET was a pair-conversion telescope. As 
such, it required a method for rejecting charged 
particles entering the detector that might other- 
wise be treated as 7-rays. This was accomplished 
by the anti-coincidence system known as the A- 
dome. The A-dome was a monolithic scintillator 
which was read out by 24 photo-multiplier tubes 
(PMTs) optically coupled to the lowest edge of the 
dome. During the EGRET beam test at the Stan- 
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ford Linear Accelerator Center (SLAC) it was dis- 
covered that the EGRET effective area {A^f / ) de- 
creased at a rate that was faster than expected at 
high energies. The decrease was determined to 
be caused by "self-veto" [2], where an otherwise 
acceptable 7-ray event is vetoed by the A-dome 
when a secondary particle associated with the elec- 
tromagnetic shower of the 7-ray in the instrument 
triggers a veto signal. The secondary particles can 
be either charged particle or X-rays which Comp- 
ton scatter within the A-dome scintillator 

Our simulation environment has allowed us to ex- 
plore the problem of self-veto with greater detail 
than was available during the calibration of the 
EGRET. Furthermore, we have simulated the in- 
strument in the environment that the astrophysics 
data were taken, specifically within close proxim- 
ity to the CGRO. Comparison between the beam 
test geometry and the flight geometry has produced 
some interesting results, primarily we find that the 
effect of self-veto is exacerbated by the EGRET 
being attached to the CGRO. 

Analysis Methods 

The simulations have a variety of parameters 
which are related to the EGRET instrument re- 
sponse. Since we are primarily concerned with the 
effect of self- veto on the A^f / we limit our discus- 
sion to the relevant parameters. Effective vetoing 
of charged particle events passing through the A- 
dome was dependent on the voltage settings of the 
readout PMTs. While in-flight, it was necessary 
to reduce the efficiency of vetoing charged parti- 
cles during calibration of the calorimeter. Thus, 
the voltage settings on the readout PMTs were ad- 
justable. To replicate this, we have implemented 
a threshold for energy measured at the PMTs for 
our simulated events. By measuring the A^jf for 
incident mono-energetic 7-rays at the same ener- 
gies measured during the EGRET beam tests we 
are able to adjust this threshold for an optimized fit 
to the EGRET reported Aeff . We have also al- 
lowed for an overall normalization factor between 
the EGRET track finding algorithm and the one 
employed in our simulations. 

The results of our optimization can be seen in Figs. 
1 and 2. This optimization was performed using 
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Figure 1: Optimization landscape with 68%, 95%, 
and 99% confidence contours. Summed veto 
threshold energy is calculated by summing the at- 
tenuated energy deposited within the A-dome. 

simulated data in the beam geometry. Every at- 
tempt was made to mimic the EGRET triggering 
algorithm and spark chamber response. However, 
as was noted by the EGRET collaboration [2], be- 
low 300 MeV there are noticeable effects on the 
track finding due to the difference in efficiencies 
between tracking layers, thus we have performed 
our fit only for energies > 300 MeV. 

For both the beam and flight geometries we gen- 
erated 800000 events at near normal incidence in 
each of the 10 standard EGRET energy bins, plus 
three extended bins (10-20 GeV, 20-50 GeV, 50- 
120 GeV), using a spectrum proportional to E^'^-^; 
this corresponds to the weighting used for genera- 
tion of the EGRET exposure maps. 

We construct scale factors for each of the 10 stan- 
dard exposure energy ranges using the following 
formula: 

(-\raccepted \ , / j^accepted \ 
i, flight i,beam \ 

y^rgenerated ] I \ i\Tgenerated j ^ 
^^i, flight / ' \^^i,beam / 

where Ni is the number of events in energy bin i, 
the subscript labels refer to the geometrical con- 
figuration of the simulation (beam test or flight), 
and the superscript labels refer to the total number 
of generated events (generated) and the events ac- 
cepted for reconstruction (accepted). These factors 
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Constructing the scaling factors as described re- 
duces the possibility of systematic effects unre- 
lated to changes in geometry. The ratio of flight 
geometry performance to beam test geometry per- 
formance should be invariant to effects unrelated to 
the change in geometry. The performance for each 
geometry has been examined in detail and there are 
no unanticipated effects that might introduce large 
systematic errors to our results. 

Results 



Figure 2: Simulated Aefj after optimizing the 
summed veto threshold energy to the beam test re- 
sults from SLAC compared to the measured values. 



are constructed assuming the same angular depen- 
dence as the published EGRET A^f / . 

Extension of the EGRET effective area beyond 10 
GeV has been done previously using the prelim- 
inary Monte Carlo mentioned earlier [8]. Our ex- 
tension uses a similar approach. The exposure map 
for the 4-10 GeV bin is used as a basis for all 
higher energy bins. The exposure maps for the 
energy bins 10-20 GeV, 20-50 GeV, and 50-120 
GeV, are generated using a scaling factor for each 
higher energy bin that accounts for the relative dif- 
ference between the number of events generated 
and accepted for reconstruction in that bin com- 
pared with the number of events for the 4-10 GeV 
bin. The following formula is used to create the 
scaling factors: 



N. 



accepted 



HE,j 



jjlight 



Tijgenera 
J^j. flight 



ted 



/A, 



accepted ^ 
9,beam 
jujgenerated 
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(2) 
(3) 



where j = 10-20, 20-50, 50-100 GeV, the super- 
scripts and subscripts have the same meaning as in 

1 J i^jaccepted . i\Tqenerated 

Eq. 1, and A^g^^jJ^ and7V|^^^„ are the num- 
ber of accepted and generated events in the beam 
configuration for the 4-10 GeV bin, respectively. 



The scale factors calculated above are applied to 
the corresponding EGRET exposure maps. All the 
scaling factors found imply a systematically lower 
Aeff at high energies. This implies that each 7- 
ray measured is more significant than previously 
thought, leading to a systematic increase in both 
the integrated flux as well as the hardness of the 
spectra measured by the EGRET. 

We have re-analyzed the EGRET data set in the in- 
ner Galaxy, 0.5° < |;| < 30° and 0.5° < |6| < 6°, 
using our re-scaled exposure maps. For compar- 
ison we have also analyzed two commonly ac- 
cepted GALPROP [9] models, 599278 [10] and 
6002029RE [11]. The GALPROP model 599278 
is constructed by assuming the observed radiation, 
gas, and cosmic -ray distributions are representa- 
tive of the Galaxy at large. This is referred to as 
the "conventional" model. 6002029RE has been 
modified to better reproduce the EGRET observa- 
tions, specifically "secondary antiproton data were 
used to fix the Galactic average proton spectrum, 
while the electron spectrum is adjusted using the 
spectrum of diffuse emission itself"[10]. This is 
known as the "optimized" model. These models 
were run in their published configurations, except 
we changed the energy and skymap binning to cor- 
respond to the EGRET skymaps. In the following, 
we denote the runs for these models as 599278EG 
and 6002029EG, respectively. Fig. 3 shows the 
E^dN^/dE^ for our 599278EG and 6002029EG 
GALPROP models as well as the re-scaled and 
original EGRET measurements. 
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Figure 3: Comparison of E^dN^/dE^ spectra for 
0.5° < |/| < 30° and 0.5° < |6| < 6°. 

Conclusions 

The re-scaled EGRET spectrum for the inner 
Galaxy shows an increased excess when compared 
to previous results. When compared with the con- 
ventional model, the EGRET excess has a reduced 

of 18.8 and 24.7 for the original and re-scaled 
EGRET measurements, respectively. For the op- 
timized model, the EGRET excess has a reduced 

of 1.9 and 3.4 for the original and re-scaled 
EGRET measurements, respectively. This analy- 
sis indicates that the GeV excess may be larger 
than previously thought. It is important to note 
that while the GeV excess is dramatic with respect 
to the conventional GALPROP model (599278EG) 
it is much reduced with respect to the optimized 
model (6002029EG). 
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